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ABSTRACT: 
Perform G floating point polynomial calculation. 
ENVIRONMENT: User Mode, AST Reentrant 
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4 .SBTTL 4 Histor 
£8: guggy SBelsinat’ Sar te eb-19 
44; 1- - Fix HL EMR. SBL 1 juty=1979 
45; 1- = Use local handler to insure that eet sour than those documented 
0 46 ; as statuses are resignalled. SBL 25-Sept-1980 


on 
.— ie | 
mm 

re 
> 
z— 
ro 
a 
—-3 


COOQooooooooooooooooooooooo 
SOoOOCCOOCSOOOOCOOOCOSOOOOOOOOoOOoOO 


COOoOoooooooooooooooOoOoooOoOooooO 


COOOCOSOCSOSOCSCOOCOCOOOC OOOO COOOOOOOOOOOOOOOOO 


oOooooooooo 
Soooooooooooeo 
Cooooooooooo 


G 1 
floating polynomial calculat 16-SEP-1984 16:21 
— vat pn Seen iobe $4 :08;81 


~SBTTL DECLARATIONS 
INCLUDE FILES: 

SCHF DEF 

SSSDEF 


; EXTERNAL DECLARATIONS: 
-EXTRN LIBSSIG_TO_RET 


> MACROS: 


+ EQUATED SYMBOLS: 


arg = 4 
degree = 8 
coeff = 12 
result = 16 


: OWN STORAGE: 


> PSECT DECLARATIONS: 


E, RD, NOWRT, LON 


C009 090909 SIN NINN NP AAA AAA AAAI 
FWN CO ODNAU SWAN (OOO NOAUES WR OWOONOULSWIN—OWo 
J 


yArLyas Macro Vv04-00 
LIBRTL.SRCJ 


-PSECT pieces Ct USR, CON, en LCL, 


LIBPOLYG.MAR; 1 


Condition handling symbols 
System symbols 


Library routine to convert 
a signal to an error return 
to caller of LIBSPOLYG. 

RO = signalled condition 


argument 

eeares of polynomial 
address of coefficient 
result of polynomial 
table 


SHR, = 
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aad -SBTTL LIBSPOLYG = Perform floating polynomial 
: FUNCTIONAL DESCRIPTION: 


LIBSPOLYG provides the functionality of the VAX hardware 
instruction POLYG to high level Language users. 


The Fhire qrerese points to a table (array) of G double 
precision floating point coefficients. The coefficient of 
the highest order term of the polynomial is pointed to 

by the table address operand, i.e. the first table element. 
The table is specified with lower order coefficients stored 
at increasing addresses. 


The evaluation is carried out by Horner's method, and the 
result is stored at the location pointed to by the fourth 
operand. The result computed is: 


OOQCoCooooooooooooo 


if d = degree 
and x = arg 
result = CLOJ+¢x*(CC1iJ+x*(Cl2J+ ... x*C£dJ)) 


The unsigned word degree operand specifies the highest 
numbered coefficient to participate in the evaluation. 
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For further detail, refer to the VAX-11 Architecture 
Handbook for the description of POLYx. 


CALLING SEQUENCE: 


status.wlc.v = LIBSPOLYG sorg.re-F. degree.rw.r, coeff.rg.ra, 
result.wg.r) 


INPUT PARAMETERS: 

arg.rg.r - argument, ‘’x'’ in polynomial 

Seitlcgte © taki of esettletente, A éndtx Giuetion 
IMPLICIT INPUTS: 

NONE 
OUTPUT PARAMETERS: 

result.wg.r - result of calculation 
IMPLICIT OUTPUTS: 

NONE 
FUNCTION VALUE: 


$S$_NORMAL 
SS$~FLTOVF 
SS$~FLTUN 
5S$~ 


successful execution 

floating overflow 

scot tae under flow ’ 
reserved operand, see VAX Architecture 
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LTUND 
ROPRAND 
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p 197 3 manual for more details | 
102 ; SIDE EFFECTS: | 
: 129 : ALL other exceptions are signalled. 
0 8 148 ; 
$00 1808 
403C 000 151 ENTRY LIBSPOLYG, “M<IV,R2,R3,R4,R5> ; Entry point, enable int. ovf. 
B08 13 3 and save R2, R3, R4, RS 
6D 16°AF 9E 000 154 MOVAB B*“HANDLER, (FP) ; Set up handler to process 
Boe 133 3; exceptions 
oc BC 08 BC 04 BC S5FD $606 157 POLYG @arg(AP), - ; perform polynomial 
QOOE 158 adegree (AP a, 3; trap on i pee to 
QOOE 159 acoef f (AP) 3; handler whic 
QOOE 160 3 unwind a return sever 
QOOE 161 3; condition in RO t 
QO0E 166 3 caller of LTB$POL YD. 
OOOE 16 
10 BC 50 7D sit: 168 MOVQ RO, @result (AP) 3; return value 
50 01 QA 0012 196 MOVZBL #1, RO 3 success status code 
04 0015 £168 RET 3 return 
0016 169 
0016 170 
0016 171 HANDLER: 
0000 0016 76 -WORD 0 
0018 17 
0018 174 i+ 
0018 175 : If the exception is one of the documented exceptions for this routine, 
Bolg 176 : call LIBSsi TO RET to return it as a status. Otherwise, resignal. 
001 177 ; Also, resignal Tf the depth is not zero. 
oo18 178 
50 O08 AC 00 oOi8 180 MOVL CHFSL_MCHARGLST(AP), RO ; Get Bochenten, vector address 
08 AO D5 Bote 181 TSTL CHFSL— _MCH_DEPTH(RO) : Is depth z 
41 12 OO1F 18 BNEQ 90$ 3 iY not. peeteaal 
51 04 aC 00 0021 18 MOVL CHESL SIGARGLST (AP), R1 ; Get signal vector address 
50 04 Ai DO 0025 184 MOVL CHFSL_SIG_NAME(R1), RO ; Get signalled condition 
048C 8F 50 B81 0029 185 CMPW = «RO, ~#SS$_FLTOVF ; Compare conditions 
gf 13 136 BEQL 10$ :; If it matches, don't resignal 
049C 8F 81 1 CMPW Re #SS$_FLTUND 
i 13 0 188 BEQL «=s«d108 
0454 8F B1 189 CMPW RO, #SS$_ROPRAND 
1¢ 13 003c 190 BEQL «=s«108 
04C4 8F 50 =—s«B1 3 191 CMP RO, #SS$_FLTUND_F 
os 12 4 135 BNEQ 5$ 
04 Al 049C BF 3C Be 19 MOVZWL #SSS_FLTUND, CHFSL_SIG_NAME(R1) ; Change fault code to trap code 
D 11 0048 19% BRBsCéd;#(S | 
0484 8F 0 ? 4 195 5$ CMP k #SS$_FLTOVF _F 
— 12 005 196 BNEQ 90% 
0s Al 048C 8F 3C 5 19 MOVZ2WL #SS$_FLTOVF CHFSL 4 51 6_MAME CRT) : 
00000000'°GF 6C FA 005A 198 10$:  CALLG (AP)> G*LIB$SIG_TOTRET ; Return signal as a status 
04 61 19 ET | 
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MOVZWL #SSS_RESIGNAL, RO 


AX/VMS Ma 
LIBRTL.SR 


it18 


04-00 
POLYG.MAR;1 


; Resignal condition 
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Page hy 


The_working set Limit was 1050 pages. 

21671 bytes (43 pages) of virtual memory were used to buffer the intermediate er 

There were 30 pages of symbol table space allocated to hold 427 non-local and 3 local symbols. 
03 source Lines were read in Pass 1, produc tng 13 object records in Pass 2. 

mac 


pages of virtual memory were used to define ros. 
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Symbol table _— ety P 90: 48:8} LIBRTL.SRCJLIBPOLYG.MAR; 1 ° (4) 
AR = (00000004 : 
CHF $L_MCHARGLST : : ; 
CHFSL_MCH_DEPTH = $ 
CHF SL_-SIGARGLST = 4 ; 
CHFSL_SIG_ = 4 ; 
COEFF = 9000000¢ ; 
DEGREE = 6 
HANDLER 9000016 R 0 ; 
LIBSPOLYG 0000000 RG ; 
LIBSSIG_10_RET eeeeeene x 0 : 
RESULT = 00000010 ; 
SS$_FLTOVF = 88 48C : 
SS$"FLTOVE _F = 4B4 ; 
SS$_FLTUN = 0000049C¢ ; 
SS$_FLTUND_F = 000004C4 : 
SS$~RESIGNAL = 00000918 : 
SS$_ROPRAND = 00000454 : 
eeeeoeeeeonoeaanaanen + P+ 
! Psect synopsis ! ; 
PSECT name Allocation PSECT No. Attributes : 
» 8 « 00000000 ( 0.) 00 ¢ OQ.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE é 
SABS$ S444 ( 0.) O01 ¢ 41.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE ; 
_LIBSCODE 000000 © 16.3 Gt 2.3 PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG ; 
$emoen enon cor een oor mawmenes ea D4 
; Performance indicators ; : 
Phase Page faults CPU Time Elapsed Time : 
Initialization 29 00:00:00.05 00:00:00.94 F 
Command processing 111 00:00:00.33 Bp a0: 05 «£8 3 
Pass 1 190 BP 282: 5-26 88289: } 36 3 
Symbol table sort 8 00: 7 0.4 00:00:01.54 F 
Pass 2 5 00:00:00.56 00:00:03.09 3 
Symbol table output 4 00:00:00.02 00:00:00.02 3 
Psect synopsis output 3 00:00:00.01 00:00:00.02 3 
Cross-reference output 00:00:00.00 00:00:00.00 3 
Assembler run totals 390 00:00:04.12 00:00:20.75 : 
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! Macro Library statistics ! 


gwar n nore neon mn ane ceeseccaa $ 


Macro- Library name Macros defined 


“$255$DUA28: CSYSLIBISTARLET. MLB;2 5 
486 GETS were required to define 5 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE =SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:LIBPOLYG/OBJ=0BJ$:LIBPOLYG MSRC$:LIBPOLYG/UPDATE=(ENH$:LIBPOLYG) 
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